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Abstract 

In the framework of the /cy-factorization approach, the production and polarization of 
prompt J/ip mesons in pp collisions at the LHC energy y/s = 7 TeV is studied. Both the 
direct production mechanism as well as feed-down contributions from Xcii Xc2 and ip' decays 
are taken into account. Our consideration is based on the color singlet model supplemented 
with the off-shell matrix elements for the corresponding partonic subprocesses. The unin- 
tegrated gluon densities in a proton are determined using the CCFM evolution equation as 
well the Kimber-Martin-Ryskin prescription. We compare our numerical predictions with 
the first experimental data taken by the CMS, ATLAS and LHCb collaborations. The esti- 
mation of polarization parameters Xg, X^ and Xg^ which determine J/ip spin density matrix 
is performed. 



PACS number(s): 12.38.-t, 13.20. Gd, 13.88. +e 



1 Introduction 

The production of charmonium states at high energies is under intense theoretical and 
experimental study p~H3]. The production mechanism involves the physics of both short and 
long distances, and so, appeals to both perturbative and nonperturbative methods of QCD. 
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This feature gives rise to two competing theoretical approaches known in the literature as 
the color singlet (CS) [I] and color octet (CO) [5] models. In the CS model, only those states 
with the same quantum numbers as the resulting charmonium contribute to the formation of 
a bound state. This is achieved by radiating a hard gluon in a perturbative process. In the 
CO model, it was suggested to add the contribution of transition mechanism from cc pairs to 
charmonium, where a charmed quark pair is produced in a color octet state and transforms 
into the final color singlet state by the help of soft gluon radiation. The CO model is based 
on the general principle of the non-relativistic QCD factorization (NRQCD) [6]. As it is well 
known, the sole leading order (LO) CS model is insufficient to describe the experimental 
data on the J / ip production at the Tevatron energies. By adding the contribution from the 
octet states and fitting the free parameters one was able to describe the data on the J/ip 
production at energies of modern colliders (see [7] and references therein). 

However, recently the next-to-leading order (NLO) [8] and dominant next-to-next-to- 
leading order (NNLO*) [9] corrections to the CS mechanism have been calculated and have 
been found to be essential in description of quarkonia production. The comparison with the 
first LHC measurements performed by the ATLAS, CMS and LHCb collaborations demon- 
strates JTU] that the NNLO* CS model correctly reproduces the transverse momentum dis- 
tributions as well as the total cross section of J/ip mesons at y/s = 7 TeV. 

The effect of high-order QCD corrections is also manifest in the polarisation predictions. 
While the charmonium produced inclusively or in association with a photon are predicted to 
be transversely polarised at LO, it has been found that their polarisation at NLO is increas- 
ingly longitudinal at high pt PELT]. Opposite, the NRQCD predicts the strong transverse 
polarization of the final state quarkonia pp. This is in disagreement with the polarisation 
measurement [12] performed by the CDF collaboration at the Tevatron, casting doubt on 
the earlier conclusion that the CO terms dominate J/ip production. 

The results of studies [8TTTT] support the predictions [T3H22] obtained in the framework of 
the fc^-factorization QCD approach [23], where investigations of heavy quarkonia production 
and polarization have own long story. Shortly, it was demonstrated [HH2T] that the experi- 
mental data on quarkonia production at HERA, RHIC and Tevatron can be well described 
within the CS model alone. The values of CO contributions obtained by fitting the Tevatron 
data appear to be substantially smaller than the ones in the NRQCD formlalism p^ | [T6 ] [2^] . 
Furthermore, the longitudinal polarization of produced J/ip mesons predicted by the kx- 
factorization is an immediate consequence of initial gluon off-shellness [H] which taken into 
account in the A/p-factorization approach^. 

In the present note we give the systematic analysis^) of first experimentl data [261428] on 
the prompt J/ip production taken by the CMS, ATLAS and LHCb collaborations at the LHC 
energy y/s = 7 TeV. Follow the guideline of previous studies [I9j|20] , in our consideration we 
will apply the CS model supplemented with the /cr-factorization approach. Two sources of 
J/ip production are taken into account: direct J/ip production and feed-down J/ip from the 
decay of other heavier prompt charmonium states like Xci, Xc2 or ip', that is in a full agreement 
with the experimental setup [26T428] . Specially we concentrate on the J/ip spin alignment 
and estimate three polarization parameters Xg, and Xg^ defining the spin density matrix 

*A detailed description and discussion of the fcT-factorization approach can be found, for example, in 
reviews [27] . 
2 See also [22|. 
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of produced J/ip mesons. As it was mentioned above, studies of polarization observables are 
useful in discriminating the CS and CO production mechanisms, 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the fc^-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 

2 Theoretical framework 

The production of prompt J/ip mesons in pp collisions at the LHC can proceed via either 
direct gluon-gluon fusion or the production of heavier P-wave states Xcj (J = 0, 1, 2) and 
S'-wave state ip', followed by their radiative decays Xcj — > J/V + 7 and if)' — > J/ip + X. In the 
CS model, the direct mechanism corresponds to the partonic subprocess g* + g* —> J/ip + g 
which includes the emission of an additional hard gluon in the final state. The production of 
P-wave quarkonia is given by g* + g* — >■ XcJ [21] and there is no emission of any additional 
gluons. The feed-down contribution from S'-wave state if)' is described by the g* + g* — > if)' ' + g 
subprocess. 

The production amplitudes of all these subprocesses can be obtained from the one for 
an unspecified cc state by the application of appropriate projection operators J(S, L) which 
guarantee the proper quantum numbers of the cc state under consideration. These operators 
for the different spin and orbital angular momentum states can be written as [I] 

J( 3 S0 = J(S = 1,L = 0) = e(S z ){p c + m c )/m l l\ (1) 

J( 3 Pj) = J(S = 1,L=1) = (p- c -m c )e(S z )(p c + m c )/m 3 / 2 , (2) 

where m is the mass of the specifically considered cc state, p c and p s are the four-momenta 
of the charmed quark and anti-quark. In accordance with the non-relativistic formalism 
of bound state formation, the charmed quark mass m c is always set equal to 1/2 of the 
quarkonium mass. States with various projections of the spin momentum onto the z axis 
are represented by the polarization vector e(S z ). 

The probability for the two quarks to form a meson depends on the bound state wave 
function ^(q). In the non-relativistic approximation, the relative momentum q of the quarks 
in the bound state is treated as a small quantity. So, we represent the quark momenta as 
follows: 

p c = p/2 + q,p 5 = p/2-q, (3) 

where p is the four-momentum of the final state quarkonium. Then, we multiply the relevant 
partonic amplitude A (depending on q) by ^(q) and perform integration with respect to q. 
The integration is performed after expanding the integrand around q = 0: 

A(q) = A\ g=0 + q a (dA/dq a )\ q=0 + (4) 

Since the expressions for A\ q= o and dA/dq a \ q= o are no longer dependent on q, they may be 
factored outside the integral sign. A term-by-term integration of this series then yields [21] 
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d>q q a *(q) = -ie a (L z )^n'(x = 0), (6) 



where lZ(x) is the radial wave function in the coordinate representation, i.e. the Fourier 
transform of ^(q)- The first term in (4) contributes only to S waves, but vanishes for P 
waves because TZp(0) = 0. On the contrary, the second term contributes only to P waves, 
but vanishes for S waves because TZ's(0) = 0. States with various projections of the orbital 
angular momentum onto the z axis are represented by the polarization vector e(L z ). The 
numerical values of the wave functions are either known from the leptonic decay widths (for 
J/ip and ip' mesons) or can be taken from potential models (for Xcj mesons). 

In our numerical calculations, the polarization vectors e(S z ) and t{L z ) are defined as 
explicit four-vectors. In the frame where the z axis is oriented along the quarkonium mo- 
mentum vector p^ 1 = (E, 0, 0, |p|), these polarization vectors read 



e*(±l) = (0,±l,i,0)/V2, e"(0) = (|p|, 0, 0, E)/m. (7) 

The states with definite S z and L z are translated into states with definite total momentum 
J and its projection J z using the Clebsch-Gordan coefficients: 

£^{J,Jz)= E (l,L z ;l,S z \J,J z )e»(S z )e u (L z ). (8) 

S Z ,L Z 

Further evaluation of all partonic amplitudes under consideration (including subsequent lep- 
tonic and/or radiative decays, of course) is straightforward and was done using the algebraic 
manipulation systems FORM [30]. We do not list here the obvious expressions because lack 
of space, but only mention several technical points. First, in according to the fc^-factorization 
prescription [23], the summation over the incoming off-shell gluon polarizations is carried 
with J2 efJ-e * u — kj.k^/ky, where k^ is the gluon transverse momentum orthogonal to the 
beam axis. In the collinear limit, when |ky| — > 0, this expression converges to the ordinary 
j2e tI e* u = —g^ u /2 after averaging on the azimuthal angle. In all other respects the evalua- 
tion follows the standard QCD Feynman rules. Second, the spin density matrix of final J ftp 
meson is determined by the momenta l\ and I2 of the decay leptons and is taken in the form 



E e'V " = 3 m + l v A - —9^ /m 2 . (9) 



This expression is equivalent to the standard one J2 e ^ e * u — —9^ u +p fl p u /m 2 but is better 
suited for studying the polarization observables because it gives access to the kinematic vari- 
ables describing the orientation of the decay plane. Third, when considering the polarization 
properties of J/ip mesons originating from radiative decays of P-wave states, we rely upon 
the dominance of electric dipole El transition^. The corresponding invariant amplitudes 
can be written as [21] 

iA(xci J/ip +i) = 9i e^Vl Xd) 4 W 4 7) . ( 10 ) 



iA{ Xc 2 -»- J/V + 7) = S2< )e f 



'111 



3 The same approach has been applied [TH] to study the T production and polarization at the Tevatron. 
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where e^ cl \ e( J /w and ejp are the polarization vectors of a corresponding spin-one particles 
and e^, c2 ) is its counterpart for a spin-two Xc2 meson, p and k are the four- momenta of the 
decaying quarkonium and the emitted photon, ^ VOL ^ is the fully antisymmetric Levita-Civita 
tensor. The dominance of electric dipole transitions for the charmonium family is supported 
by the experimental data taken by the E835 Collaboration at the Tevatron [32]. Since the 
electromagnetic branching ratio for Xco J/ip + j decay is more than an order of magnitude 
smaller than those for Xd an d Xc2, we neglect its contribution to J /if) production. As the 
if)' — )■ J /if) + X decay matrix elements are unknown, these events were generated according 
to the phase space. 

The cross section of J /if) production at high energies in the fc^-factorization approach is 
calculated as a convolution of the off-shell partonic cross section and the unintegrated gluon 
distributions in a proton. The contribution from the direct production mechanism can be 
presented in the following form: 



a(pp -> J/if) + X) = f 1 -f g (x 1 ,k 2 1T ^ 2 )f g (x 2 ,k 2 2T ^ 2 ) 

J lb7T{XiX 2 S) z 

x \M(g* + g* ->■ J /if) + g)\ 2 dp^dk 2 1T dkl T dydy, 



2tt 2tt ' 

where f g (x, k^, /i 2 ) is the unintegrated gluon density, pr and y are the transverse momentum 
and rapidity of produced J /if) meson, y g is the rapidity of outgoing gluon and s is the pp 
center-of-mass energy. The initial off-shell gluons have a fraction x\ and X2 of the parent 
protons longitudinal momenta, non-zero transverse momenta k\T and k2T (kf T = — k\ T ^ 0, 
^2T = ~^2T 7^ 0) an d azimuthal angles <pi and <f> 2 . For the production of Xcj mesons via 
2—7-1 subprocess above we have 

r 2tt 

a(pp -»> XcJ + X) = / -/ 9 (x 1 ,k 2 T ,/i 2 )/ 9 (x 2 ,k 2 T ,/i 2 )x 

J X1X2S 1 (13) 

x \M(g* + g* -+ Xc j)\ 2 dkj T dkl T dy 

ATI 271 

where T is the off-shell gluon flux factor. In the present analysis we set it to be equal to 
T = 2s, where s is the energy of partonic subprocess. In (12) and (13), \M.(g* + g* — > 
J/ip + g)\ 2 and \M.{g* + g* — > Xcj)\ 2 are the corresponding off-shell matrix elements squared 
and averaged over initial gluon polarizations and colors. The production scheme of if)' meson 
is identical to that of J /if), and only the numerical value of the wave function |7£(0)| 2 is 
different (see below). 

In the numerical calculations we have tested a few different sets of unintegrated gluon 
distributions involved in (12) and (13). First of them (CCFM set AO) has been obtained [33] 
from the CCFM equation where all input parameters have been fitted to describe the proton 
structure function F 2 (x,Q 2 ). Equally good fit of the F 2 data was obtained using different 
values for the soft cut and a different value for the width of the intrinsic ky distribution 
(CCFM set BO). Also we will use the unintegrated gluons taken in the Kimber-Martin- 
Ryskin (KMR) form [3l]. The KMR approach is a formalism to construct the unintegrated 



parton distributions from well-known conventional ones. For the input, we have used recent 
leading-order Martin-Stirling-Thorn- Watt (MSTW) set [35]. 
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The multidimensional integrations in (12) and (13) have been performed by the means of 
Monte Carlo technique, using the routine VEGAS |36j. The full C++ code is available from 
the author on request^. 

3 Numerical results 

We now are in a position to present our numerical results. First we describe our input 
and the kinematic conditions. After we fixed the unintegrated gluon distributions, the cross 
sections (12) and (13) depend on the renormalization and factorization scales \Lr and fip. 
Numerically, we set = m 2 + and fi 2 F = s + Q^, where Qt is the transverse momentum 
of initial off-shell gluon pair. Note that the choice of hr is the standard one for studying of 
the J/if> production whereas the special choice of /xp is connected with the CCFM evolution 
(see |33j). Following to [37J, we set mj/^ = 3.097 GeV, m Xcl = 3.511 GeV, m Xc2 = 3.556 GeV, 
m^i = 3.686 GeV and use the LO formula for the coupling constant a s (/i 2 ) with rif = 4 quark 
flavours at Aqcd — 200 MeV, such that a s (M|) = 0.1232. The charmonia wave functions 
at the origin of coordinate space are taken to be equal to \TZj/^(0)\ 2 = 0.0876 GeV 3 [37], 
\TZ' x (0)\ 2 = 0.075 GeV 5 [38], \TZ^(0)\ 2 = 0.0391 GeV 3 [37]. According to [37], the following 
branching fractions are used: B(xci ~^ J + l) = 0.356, B(xc2 ~> J/i> + 7) = 0.202, 
B(ip' -+ J /if) + X) = 0.561 and B(J/ip -+ /i+ii~) = 0.0593. 

The results of our calculations are presented in Figs. 1 — 3 in comparison with the CMS, 
ATLAS and LHCb data [2"6T - T2"8"] . The solid, dashed and dash-dotted curves correspond to 
the results obtained using the CCFM AO, B0 and KMR gluon densities, respectively. Ev- 
erywhere, we separately show the contribution from the direct production mechanism taken 
solely (dotted curves). In this case we apply the CCFM AO gluon density for illustration. 
It is clear that sole direct production is not sufficient to describe the LHC data. However, 
we obtain a good overall agreement of our predictions and the data when summing up the 
direct and feed-down contributions. The latter is important and production of J /if) mesons 
via radiative decays of Xcj and if)' mesons even dominates over the direct contribution at 
large transverse momenta. The reason can be seen in the fact that the production of Xcj 
states refers to much lower values of the final state invariant mass and therefore effectively 
probes small x region, where the gluon distributions are growing up. The dependence of 
our numerical results on the unintegrated PDFs is rather weak and the CCFM and KMR 
predictions are practically coincide. The difference between them can be observed at small 
Pt or at large rapidities probed at the LHCb measurements. 

Computations [7] performed in the framework of NRQCD, where CO contributions are 
taken into account, can also explain at satisfactory level the shape and the absolute normal- 
ization of the measured J /if) cross- sect ions. However, as it was mentioned above, they predict 
a substantial transverse component for the polarisation of J /if) mesons at large pt which is 
not supported by measurements. We find that in the framework of the /cT-factorization ap- 
proach no need for a CO contributions in the description of J /if) production at the LHC. 
From the other side, the account of high-order corrections to the CS cross sections calculated 
in the collinear QCD factorization also leads to the significant improvements in description 
of the data: the upper bound of the NNLO* CS predictions is very close [TT] to the mea- 

4 lipatov@theory.sinp. msu.ru 
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Source 


Xe (HX) 


A^ (HX) 


\o<j> (HX) 


Xe (CS) 


A<£ (CS) 


Ae»0 (CS) 


Direct 


-0.15 


-0.09 


0.01 


0.20 


-0.22 


-0.01 


r ecu. LIU W 11 


n i q 

u. ly 


n i a 

U. 11 


n nn 


n ^ 

u.oo 


n dq 

u.uj 


n nn 

u.uu 


Total 


-0.07 


-0.03 


0.01 


0.24 


-0.14 


-0.01 



Table 1: The polarization parameters of prompt J /if) mesons calculated in the kinematical 
region of CMS and ATLAS measurements [26,27J. The CCFM AO gluon density is used. 

surements [261428] and agree much better (compared to the LO CS results) with the kx- 
factorization calculations which incorporates a large part of collinear high-order corrections 
at LO level. 

Note that the calculated cross sections of feed-down contributions from the P-wave states 
are free from singularities at small transverse momenta. This contrasts with the collinear 
QCD factorization predictions, which are either unphysical or even divergent. 

Now we turn to the the J /if) polarization. In general, the spin density matrix of a vector 
particle depends on three parameters Xg, X^ and Xg^ which can be measured experimentally. 
So, the double differential angular distribution of the J /if) — > ji + decay products reads [39] 

- — ° ,, ~ 1 + X e cos 2 9* + Xs sin 2 9* cos 20* + Xgs sin 29* cos 0*, (14) 
a cos v*d(f>* 

where 9* and 0* are the polar and azimuthal angles of the decay lepton measured in the 
J /if) rest frame. Since the polarization parameters Xg, X^ and Xg^ (which greatly affects on 
the cross sections) are not determined yet at the LHC, the results of measurements per- 
formed by the CMS, ATLAS and LHCb collaborations have been presented in a different 
ways. So, in the ATLAS analysis [27] the unknown J / if) polarization has been treated as an 
additional source of systematic uncertainties. Contrary, the CMS and LHCb collaborations 
quote their measurements [261128] for different polarization scenarios: unpolarized (Xg = 0), 
full longitudinal polarization (Xg = —1) and full transverse J /if) polarization (Xg = 1) in the 
Collins-Soper or the helicity framed. Below we estimate the polarization parameters Xg, X^ 
and Xg<f) in a whole kinematical regions regarding the CMS, ATLAS and LHCb measure- 
ments. Our evaluation is generally followed the experimental procedure. We have collected 
the simulated events in the specified bins of J /if) transverse momentum and rapidity y, 
generated the decay lepton angular distributions according to the production and decay ma- 
trix elements, and then applied a three-parametric fit based on (14). The estimated values of 
polarization parameters Xg, X^ and Xg^ in the helicity (HX) and Collins-Soper (CS) frames 

5 The experimental data points in Figs. 1 and 3 correspond to the unpolarized scenario. 
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Source 


Xe (HX) 


A^ (HX) 


\e<t> (HX) 


Xe (CS) 


A^ (CS) 


Xq^ (CS) 


Direct 


-0.03 


-0.13 


0.17 


0.19 


-0.22 


-0.03 


r ecu. U.U W 11 




nil 

U.ll 


n i ^ 

u. io 




u.uo 


u.uo 


Total 


0.03 


-0.07 


0.16 


0.26 


-0.14 


-0.01 



Table 2: The polarization parameters of prompt J/tp mesons calculated in the kinematical 
region of LHCb measurements [28J. The CCFM AO gluon density is used. 

are listed in Tables 1 and 2. We find that these parameters are the same in the kinematical 
regions covered by the CMS and ATLAS collaborations. In order to study the production 
dynamics in more detail, we separately show contributions from the direct and feed-down 
mechanisms. The latter, of course, change the polarization of final J /if} mesons predicted 
by the direct production mechanism [Tj5] but this effect is not well prononced due to overall 
integration over J /if) transverse momentum. Note that the qualitative predictions for the 
J /if) polarization are stable with respect to variations in the model parameters. In fact, 
there is no dependence on the strong coupling constant and unintegrated gluon densities, 
i.e. two of an important sources of theoretical uncertainties cancels out. Therefore future 
precise measurements of the polarization parameters at the LHC will play crucial role in 
discriminating the different theoretical approaches. 

4 Conclusions 

We have investigated prompt J / if) production in pp collisions at the LHC energy y/s = 
7 TeV within the framework of the /^-factorization approach. Both the direct production 
mechanism as well as feed-down contributions from Xci, Xc2 and if)' decays are taken into 
account. Our consideration is based on the color singlet model supplemented with the off- 
shell matrix elements for the corresponding partonic subprocesses. The unintegrated gluon 
densities in a proton are determined using the CCFM evolution equation as well the Kimber- 
Martin-Ryskin prescription. We have obtained well agreement of our calculations and the 
first experimental data taken by the CMS and ATLAS collaborations when summing up the 
direct and feed-down contributions. The dependence of our predictions on the unintegrated 
gluon densities appears at small transverse momenta and at large rapidities covered by the 
LHCb experiment. We have demonstrated also that in the framework of the /^-factorization 
there is no room for a color octet contributions for charmonium production at the LHC. 

The estimation of the polarization parameters X$, X^ and Xe<f> which determine the J/if> 
spin density matrix is given. The future experimental analysis of the quarkonium polarization 
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at the LHC turned out to be very important and informative for discriminating the different 
theoretical models. 
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Figure 1: The double differential cross sections da/dydpr of prompt J /if) production at 
y/s = 7 TeV compared to the CMS data [26j. The solid, dashed and dash-dotted curves 
correspond to the results obtained using the CCFM AO, CCFM BO and KMR gluon densi- 
ties, respectively. The dotted curves represent the contribution from sole direct production 
mechanism calculated with the CCFM AO gluon distribution. 
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Figure 2: The double differential cross sections da/dydpr of prompt J/ip production at 
y/s = 7 TeV compared to the ATLAS data [27] . Notation of all histograms is the same as 
in Fig. 1. 
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Figure 3: The double differential cross sections da/dydpT of prompt J/ip production at 
a/s = 7 TeV compared to the LHCb data [28J. Notation of all histograms is the same as in 
Fig. 1. 
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